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ABSTRACT 

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out 
deep extragalactic surveys, one of whose aims is to establish spectral energy distributions 
(SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength re- 
gion. We report observations of the (IR/SMM) emission from the Lockman North field (LN) 
and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy im- 
ages in the wavelength range covered by Herschel generally represent a blend with contribu- 
tions from neighboring galaxies, we present sets of galaxies in each field especially free of 
blending at 250, 350, and 500 fim. We identify the cumulative emission of these galaxies and 
the fraction of the far infrared cosmic background radiation they contribute. Our surveys re- 
veal a number of highly luminous galaxies at redshift z < 3 and a novel relationship between 
infrared and visible emission that shows a dependence on luminosity and redshift. 

Key words: IR galaxies: spectral energy distributions: galaxy luminosities 



1 INTRODUCTION 



The Herschel Space Observatory (Pilbratt et al. 2010) has 
opened wide astronomical access to the far-infrared/submillimeter 
(FIR/SMM) spectral range. With the Spectral and Photometric 
Imaging Receiver (SPIRE) (Gri ffin et al.|2010) , deep cosmologi- 
cal surveys are studying galaxies out to redshifts of order z ~ 3, 
reaching back to epochs when the Universe was only a few billion 
years old. 

A primary motivation for these surveys, as well as those un- 
dertaken with the Photodetector Array Camera & Spectrometer 
(PACS) aboard Herschel, (Poglitsch et al. 2010 1 is to gain improved 
spectral energy distributions (SEDs) of astronomical sources. With 
an appropriate redshift, integration of the flux densities demarcated 
by the SEDs permits derivation of rest-frame luminosities, star- 
formation rates, and other physical properties of galaxies. 

Among the first observations undertaken by SPIRE in the Her- 



schel Multi-tiered Extragalactic Survey (HerMES ) project ( Oliver 



etal. 20101 have been surveys of galaxies in GOODS-N and North- 
ern portions of the Lockman Hole (LN) field [see Oliv er et al.| 
(2010a) for a description of these early observations.] Source con- 
fusion, as defined and discussed in detail by Take uchi, T. T. &] 
|Ishii, T. T]j2004| , results in blending of far-infrared sources and 
complicates the analysis of survey data. In light of the large de- 
gree of source blending expected at SPIRE wavelengths, novel op- 
tions for source extraction have been pursued [(Roseboo m et al.] 
|20"T0l >, (Smith et al.l[20To) , ( |Bethermin et al.|2010| >]. Rather than 
looking for sources based on SPIRE intensity maps alone or rely- 
ing on traditional source detection and extraction techniques for the 
SPIRE data, which are heavily affected by confusion, Rose boom et| 
[all pOTol ) measure the SPIRE flux at the position of known 24 fim 
sources using a linear inversion technique to account for source 
blending. The rationale for this is provided by the results of the 
Balloon-borne Large Aperture Submillimetre Telescope (BLAST) 
extragalactic survey ( jMarsden et al. 2009 1, which showed that the 
24 yum and the FIR flux densities are at least statistically correlated. 



2 PRIMARY AIMS 

The aims of this paper are twofold; our primary aim is to de- 
rive spectral energy distributions for distant galaxies observed by 
SPIRE. Before this can be achieved, however, a robust way of iden- 
tifying sources least affected by confusion and blending must be 
devised. 

The GOODS-N catalogue of |Roseboom et al.|p0T0) provides 



a cross-identification (XID) of FIR/SMM flux density at 250, 350 
and 500 fim with 1951 possible 24 jim counterparts having mini- 
mum flux densities of 20 yuly. Many of the identified 24 jim galaxies 
are further cross-identified with ultraviolet, optical, near-infrared 
(NIR) and radio counterparts. The survey covered a 12.3 x 18.6 arc 
minute strip, corresponding to ~230 arcmin 2 

The SPIRE beam diameters at full-width-half maximum 
(FWHM) respectively measure 18.1, 25.2 and 36.9 arcsec at 250, 
350 and 500 fim. For present purposes, we take the beams to be 
close to circular; their ellipticity varies from pixel to pixel, but is 
approximately 1.08 ± 0.05, the longer direction lying in the space- 
craft horizontal direction, parallel to the ecliptic plane ( B ernhard| 
Schulz private communication I. The beam at 500 fim thus has an 
area ~0.3 arcmin 2 . With 1951 possible 24 yum sources, we can ex- 
pect a typical 500 fim beam to contain 2.5 possible sources. At 250 
fim the crowding is a factor of 4 less severe, but still appreciable. 
An example of the crowded source distribution is seen in Figure 
[T] where the SPIRE beam outlines are overlaid on a patch of the 
GOODS-N field centered on a 24 jim source. 

In the Northern Lockman region, a 40.1' x 36.2' segment 
of the sky yielded 6316 possible 24 jim counterparts with mini- 
mum flux densities of 50 /jly in an area subtending ~1500 arcmin 2 
— again corresponding to more than 1 potential 500 jim source 
per beam. Identification of sources least affected by confusion 
and blending is therefore important. These sources have also been 
cross-identified at multiple wavelengths and assigned photometric 
redshifts as detailed in Stra zzullo et al.|p010| >. 



3 DATA 

Using the existing high spatial resolution Spitzer 24 fim data and 
the known far-infrared instrumental point response function (PRF) 
as inputs, Roseb oom et al.| p010) determined best-fit 250, 350, and 
500 fim fluxes by a procedure they detail in their paper. In each 
SPIRE wavelength band, their tabulated cross-identifications pro- 
vide both their best estimate of the flux density F v and the flux 
density PRF Y in a PRF-convolved map centered on the position of 
an associated 24 yum source. They make no assumptions about a 



1 Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 



ipation from NASA. 
2 hermes.sussex.ac.uk 
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proportionality between 24 pm and FIR flux densities, but assume 
that SPIRE sources will only be detected at positions of 24 pm 
sources. For each source, ancillary data at other wavelengths are 
included, as well as several flags that can be used to identify de- 
generate cases. Of particular importance is the redshift of the asso- 
ciated 24 fim source, which enables derivation of rest-frame SEDs 
and thus source luminosities. Wherever we refer to F v and PRF V in 
the remainder of this paper these quantities are to be thought of as 
those defined by Rose boom et al.| ( [2010| >. 

The flux density in a PRF-convolved region centered on the 
position of the 24 pm galaxy associated with each far-infrared 
source represents the system response not only to the flux density 
attributed to this source (referred to as F v \^) but also to contribu- 
tions from nearby sources. It thus provides a measure of blending 
characterizing each source. We define a purity index IT,! for each 
source as the ratio 



n A = fjprfJ 



(i) 



where A specifies the wavelength band, 250, 350, or 500 pm, and 
PRF V \ , as supplied in the XID catalogs by Roseboom et al. 2010 
is the PRF-smoothed flux density at the position of the source. 3 A 
high value of Tl A indicates low confusion and blending in wave- 
length band A; a low value indicates high blending. The fractional 
contributions by ambient sources to the PRF-convolved flux den- 
sity within a PRF is simply (1 - II,!). In principle, the purity index 
must assume values < T1 A < 1. In practice, however, the value 
F y has been calculated on top of a. locally determined background, 
whereas the PRF V values do not take a local background variation 
into account. This can result in Tl A > 1 when a negative local back- 
ground or "baseline" has been used. 



4 THE STATISTICS OF PURITY INDICES IN GOODS-N 
AND LOCKMAN NORTH 

We have found it useful to identify sources whose purity indices 
respectively exceed Tl A measures of 0.7, 0.5, and 0.3 at 250, 350, 
and 500 pm. We say that a source is secure at each wavelength if it 
meets this criterion. If the source is found to be secure in all three 
wavelength bands, we call it triply secure. High-purity sources are 
of special interest because their isolated nature makes them less 
susceptible to blending by neighbors. It is these sources which will 
be especially useful for follow up studies with other instruments 
and also may provide confirmation of the deblending approach 
used. It should be noted, however, that a flux density estimate from 
a highly pure source might still be inaccurate if there is signifi- 
cant contribution from an infrared source that is not observed at 24 
pm. Furthermore, sources with low purity do have well-defined de- 
blending solutions and hence well-characterized flux densities and 
uncertainties from Roseboom et al. (2010). In crowded fields, the 
true flux density is very likely described by this characterization, 
although the margins of uncertainty tend to be large. Our choice 
of purity criteria is somewhat subjective but offers a reasonable 
compromise for extracting relatively reliable SEDs despite source 
blending. Understandably, these criteria may be expected to vary 
depending on the type of information an astronomer expects to ex- 
tract from the survey data. 

3 Our PRF,, |,i corresponds to the quantity d in equation (2) of Roseboom 
et al. 2010, convolved with the point response function centered on the pri- 
mary source whose flux is F r \ A . The entries in the XID tables list our F v \ A 
as F(A), and our PRF V \ A flux density as PRF(A). 



When a reliable redshift is available, the single most important 
quantity that can be determined from an SED is source luminos- 
ity. With this, one can begin discussing the luminosity distribution 
at specific redshifts, as well as luminosity evolution as a function 
of redshift, particularly among ultraluminous galaxies that emit the 
dominant fraction of their energy in the infrared. However, to obtain 
a reliable SED and thus a reliable luminosity, we require sources 
whose flux densities are well determined at all three SPIRE wave- 
lengths in order to optimally constrain the flux density defining the 
broad wavelength region around peak emission. 

To explain the consequences of our choice of purity criteria in 
this context, we may consider a toy model which, as pointed out in 
Section 1, will respectively exhibit an average number of sources 
n A ~ 2.5, 1.225, and 0.625 per GOODS-N beam, at 500, 350, and 
250 pm. Let us inject an additional source into such a beam and call 
it the primary source. If all the sources involved are equally bright, 
on average, i.e., make equal contributions to the PRF-smoothed flux 
density, the purity of the primary source will be IL, = (1 + n A )~ l , 
i.e. 0.29, 0.45, and 0.62, respectively at 500, 350 and 250 pm. Half 
the sources in each waveband will have purities higher than these 
purity cuts, and half lower. 

Turning now to our preferred adoption of purity cuts of 0.3, 0.5 
and 0.7 at 500, 350, and 250 pm, we see that they assure two prop- 
erties: (i) that they yield sources whose purities are above average at 
all three wavelengths, and (ii) that the fraction of sources with pu- 
rity above the cut is roughly comparable at all three wavelengths — 
a balance, which is important to assure a well-defined SED. Table 
1, described below, confirms these traits for the GOODS-N sample. 
It shows that a fraction f A = 0.23 of the sources has purity exceed- 
ing 0.7 at 250 pm, a fraction 0.32 exceeding purity 0.5 at 350 pm, 
and a fraction 0.36 exceeding purity 0.3 at 500 pm. These fractions 
cluster around a value of 0.3, thus lending roughly equal weight to 
the flux density in each waveband in the determination of the SED. 

In Table 1, we list the fraction of sources in GOODS-N and 
LN whose purity indices lie above certain cuts. We permit these 
indices to slightly exceed a value of one, with a cut-off of n = 1.1 
in GOODS-N and n = 1.2 in LN. These relaxed upper limits are 
designed to allow the inclusion of detections with a significant local 
background which otherwise have all the earmarks of being secure. 

In the GOODS-N region the XID catalog lists 183 sources 
observed at all three wavelengths and with known redshift, 16 of 
which are triply secure. If we remove the upper limit n = 1.1, the 
number of triply secure sources rises to 59. In LN there are 633 
sources with detections at all five wavelengths with known red- 
shift, 165 of which are triply secure; this number increases to 287 
if we remove the upper limits on purity. Although the numbers of 
these galaxies are quite modest, they nevertheless yield informa- 
tive statistics on the luminosities and luminosity distributions of 
galaxies observed out to redshifts z ~ 3. These will be discussed in 
Section 6. 

The larger aperture of the Herschel telescope and the higher 
spatial resolution this enables have permitted the SPIRE surveys to 
reach depths beyond those attained by BLAST. Nevertheless, Mars- 
den et al. (2009 ) succeeded in acquiring reliable measurements of 
stacked source flux densities at comparable wavelengths. Their re- 
sults indicate surface brightnesses of 8.60 ± 0.59, 4.93 ± 0.34, and 
2.27 ± 0.20 nW mT 2 sr 1 at 250, 350, and 500 pm respectively. 
These stacked source flux densities represent the major compo- 
nent of the cosmic infrared background (CIB) measured by the 
Cosmic Background Explorer's Far-Infrared Absolute Spectrom- 
eter (FIRAS) to be 10.4 ± 2.3, 5.4 ± 1.6, and 2.4 ± 0.6 nW m~ 2 
sr 1 at 250, 350, and 500 pm respectively ( jFixsen et al.|l998"| l. To 
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Table 1. Fraction of detected SPIRE sources in GOODS-N and LN with a 
IL value within the range indicated in the top row. The columns marked 
"detections" denote the total number of SPIRE sources detected at a given 
wavelength in the current HerMES survey. 
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investigate the extent to which this background is resolved with 
Herschel, we summed the estimated flux densities for our secure 
SPIRE sources in the deepest field (GOODS-N) and then attributed 
this flux density to the entire survey region of 230 arcmin 2 . Using 
the flux densities for sources contained in the XID catalog o f|Rose-| 
|boom et al.|p0T0) , our cumulative surface brightness for GOODS- 
N at 250, 350, and 500 fim is 1.49, 0.70, and 0.41 nW nT 2 sr _1 
or 14%, 13%, and 17% of the estimated CIB. At 250 and 350 yum 
these values are within 1 cr of those corrected for blending and in- 
completeness by Oliver et al. (2010al. 



5 SPECTRAL ENERGY DISTRIBUTION OF THE 
SECURE SOURCES 

Figures [2] and [5] exhibit SEDs for triply secure sources in LN and 
GOODS-N. For LN sources, we have set an additional criterion 
for inclusion, namely that they have observed flux densities also at 
PACS wavelengths of 100 and 170 jim. Along with examining the 
observed SED, we show a fit using starburst models developed by 
Siebenmorgen & Kriigel (20071 (S&K). These models are based 
on a nuclear concentration of massive young stars embedded in a 
matrix of gas and dust referred to as "hot spots". S&K use a five 
parameter SED fit which incorporates a variable nuclear bulge size 
with old and new stellar components as well as the effects of dust. 
They provide their models in the form of a library of 7000 SEDs 
available as text files online. 4 By using their models, we are able to 
not only find realistic intrinsic luminosities, but also star formation 
rates (SFR) for highly luminous sources at high redshift for which 
the Kennicut infrared - SFR relations apply (Kennicutt 1998 1. 

The S&K model fits observations quite well, although some 
of our SEDs exhibit considerable deviations from the data at vis- 
ible and near infrared wavelengths. This is largely due to variable 
shielding of starlight by dust, see section 6 below. 

In Figures [2] and [3] we have focused on high redshift (z>0.5) 
sources as these are of greatest interest to star formation history. 
Tables [2] and [3] show the luminosities and star-formation rates of 
galaxies whose SEDs are shown in Figures 2 and 3. The models 
appear capable of correctly representing the galaxies we observe, 
most of which are highly luminous and likely to be starbursts. Like 




189.1 2 IS&rlB" 189.18 

RA 



Figure 1. A GOODS-N SPIRE and 24 /im source with equatorial coordi- 
nates (J2000.0) 189.1506 / 62.2824 shown as a diamond with nearby 24 /an 
sources. Known ambient 24 /im sources are shown as squares with nested 
symbols. Plus signs indicate sources with known redshifts and triangles in- 
dicate sources with unknown redshift. The spatial distribution of 24 fim 
sources is based on the XID information in Roseboom et al. ( 20101. Solid 
circles represent the full-width-half-maxima (FWHM) of the Airy profile 
PRF at the three SPIRE wavelengths and dashed circles represent the sec- 
ond Airy minimum. Coordinates are in J2000. 

most current models, however, vastly different dust masses and nu- 
clear sizes are able to yield similar SEDs, as witnessed by the large 
differences in dust masses assigned to some galaxies listed in Ta- 
bles]^ and [3] that have nearly identical redshifts and luminosities. 
This is not surprising because the models are only required to pro- 
vide sufficient dust to convert most of the visible and ultraviolet 
radiation produced in the starburst into infrared radiation at the ob- 
served temperature. If this criterion is satisfied, the models produce 
roughly correct SEDs. 

While these figures and tables emphasize the relatively few 
triply-secure SEDs among high redshift objects, it is important to 
note that for many statistical trends, triple-security is not neces- 
sary. In Figure [4] we show the luminosities we have determined as 
a function of redshift, based on the S&K models for all sources in 
GOODS-N and LN that have detections of any kind (secure or not) 
at all three SPIRE bands. While some of the flux densities, espe- 
cially at 500 /im, are quite uncertain (see Figures|2]and[3](, this will 
not greatly affect the overall luminosity distribution. We estimate 
the source luminosity uncertainty to be ~20%. The minimum de- 
tectable luminosities follow the expected trend of increasing with 
the square of luminosity distance (which we determined based on 
a A = 0.7, H = 70 km s~' Mpc~' cosmology). 

Throughout this paper we refer to multiple luminosities, 
namely total bolometric luminosity and infrared luminosity inte- 
grated over 8 - 1000 pm [as used in star formation estimates by 
|Kennicutt| ( fl998] >]. For our purposes, the difference between the 
two is small as the majority of a luminous star-forming galaxy's 
energy is emitted in the infrared. In galaxies with L\ R > 10 H 5 L Q , 
|Buat et al.| l |2010) find that ~95% of the total star formation rate is 
accounted for by L IR . Nonetheless, we explicitly differentiate be- 
tween the two when the distinction is significant. 



6 THE NATURE OF THE ULTRALUMINOUS SOURCES 



4 The S&K SED library of models is available at 
http://www.eso.org/~rsiebenm/sb_models/ 



As Figure 4 attests, our surveys reveal a number of highly luminous 
sources, mainly at redshifts between z = 2.5 and 3.2. Discussing 
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all of these in detail is beyond the scope of the present paper, if for 
no other reason than that information on these sources is still quite 
modest. Nevertheless, we here discuss four of these sources to show 
that their luminosities appear to be intrinsic. Neither gravitational 
lensing nor blending from neighboring sources appear to contribute 
significantly to the observed luminosities. This may not be true of 
all of the luminous SPIRE sources, but it appears to be so for at 
least those about which we have the most information right now. A 
number of other uncertainties also warrant comment. 

(i) For most of the sources observed in LN only photometric 
redshifts are available. For the more luminous sources shown, es- 
timated redshift errors range from ~ 6 to 17% in the catalogs of 
|Roseboom et al.H201(ty . Because calculated luminosities are pro- 
portional to (1 + z) 2 , a 17% error at z = 2.5 can lead to a luminosity 
error in the range of ~ 23 - 26%. For sources with established 
spectroscopic redshifts, the uncertainty is far smaller and can be 
neglected. 

(ii) A second area of concern is discrimination between intrin- 
sically luminous sources and sources that may be lensed by fore- 
ground galaxies to merely appear luminous. For some purposes, 
e.g. accounting for sources contributing to the diffuse background, 
this distinction may not be of primary importance; but for charting 
the luminosities of distant sources their intrinsic luminosities need 
to be determined. 

(iii) Because our deep surveys are highly sensitive, they re- 
veal a large number of faint sources. This leads to potential mis- 
identification of sources. It can also lead to blending. But seven of 
the eleven ultraluminous sources in Figure 5 and Table 4 are triply 
secure; the remaining four are doubly secure within observational 
uncertainties. Severe blending is thus unlikely. 

(iv) For some of the sources, the available data points strad- 
dle but do not directly constrain the region where the SED reaches 
a maximum (see Figure 5), so that our SED-derived luminosities 
could be over-estimated. 

At the present stage of data reduction we are not yet in a posi- 
tion to account for all these uncertainties. However, for the cited 
sources in GOODS-N and for a few of the sources in LN reli- 
able spectroscopic redshifts are available. Among four of these we 
were able, below, to search for potential lensing, assess a degree of 
blending, and justify our confidence in their derived luminosities in 
some detail. These sources are referred to by letters corresponding 
to their designations in Table|4] Note that in Table|4]we give bolo- 
metric luminosities whereas below we quote infrared luminosities. 

The LN source (j) at equatorial position (J2000.0) 161.554052 
/ 58.788592, is characterized in the Sloan Digital Sky Survey 
(SDSS) as a well-isolated circular source identified as a quasi- 
stellar object (QSO). The visible spectrum leaves no doubt about 
the redshift z = 3.037 determined by the strong Ly-a line. The visi- 
ble continuum flux in this spectrum is ~ 5xl0~ 17 ergcirr 2 s _I A -1 , 
equivalent to ~ 4 x 10~ 5 Jy, stretching down to ~ 1,000 A, for a 
total flux of roughly 1.2x 10~ 15 W nT 2 , which is comparable to the 
infrared flux observed from the source. The IPAC extragalactic data 
base lists the optical source as brighter, by more than a magnitude, 
tha n any other source within a radius of an arc minute. |Trouille et] 
aT1 j2008f list a 0.5-2 keV X-ray flux of 38.4 x lO" 18 W nr 2 , and a 



2-8keV flux of 22.4 x 10~ I8 W irT 2 for the source, jointly about 20 
times lower than the visible flux. Owen & Morrison ( 2008 1 detected 
a 20-cm continuum radio flux density of 64 yidy from this source, 
within an apparent size < 1.5 arcsec FWHM, ruling out obvious 
ambient emission that might have indicated lensing. The XID cata- 
logs of Roseboom et al. (2010} show a relatively weak neighboring 
SPIRE source at a distance of 26 arcsec, and a source comparable 



in infrared brightness to LN 4241 but at a separation of ~ 34 arcsec. 
Because of their relatively large displacement, these and other am- 
bient galaxies are unlikely to contribute appreciably to the SPIRE 
flux densities assigned to LN (j) by Roseboom et al. (2010). All this 
gives confidence that both the visible and the infrared flux come 
from the same source, that there is no lens magnification, and that 
we are indeed dealing with an un-lensed ultraluminous source with 
infrared luminosity ~ 1.8 x 10 13 L o . The simplest explanation for 
these data is that we are viewing a QSO with a surrounding torus 
along a sightline coinciding with the torus axis. The visible light 
reaches us directly along this axis; the infrared emission comes 
from dust heated in part by the QSO and possibly also by massive 
star formation. 

The GOODS-N source (c), with equatorial coordinates 
(12000.0) 188.990097 / 62.17342, is cited by |Barger et al~|p008l l 
as having spectroscopic redshift z = 3.075. The 24 fim flux den- 
sity listed in the Barger catalog is 109 jiiy. The rest-frame 2- 
8keV luminosity assigned to this source by |Trouille et~aL]p008] > 
is 3.464 x 10 36 W ~ 9.01 x 10 9 L o , but their paper provides no rest- 
frame 0.5-2keV luminosity. The nearest SPIRE source listed in the 
XID catalogs lies at a distance of 35 arcsec, where its contribution 
to our primary galaxy's flux listed in the XID catalog can at most 
be minor. Our primary source also displays weak X-ray fluxes, but 
|Morrison et al.| j2010) list no 20-cm source within 3.5 arcmin. In- 
tegrating the flux densities indicated by the fitted SED in Figure|5] 
leads to an infrared luminosity of 6.0 x 10 I2 L o . 

One of our ultraluminous sources [LN (f)] has previously been 
discussed by Po lletta et al.| ( |2006| l. They observed the Lockman 
SWIRE source at (J2000.0) 161.041521 / 58.87355 with Chandra 
and detected a flux of 2.7 ± 1.1 x 10~ 15 erg cirT 2 s -1 in the 0.3 
- 8 keV range. The Spitzer 24 fjm flux is 4.0 mjy, strong for a 
source at spectroscopic redshift 2.54, and much brighter than any- 
thing within an arcminute of its location. Because of its initial de- 
tection by Spitzer, the authors characterize the source as an infrared 
selected Compton-thick AGN on the basis of the rest-frame hydro- 
gen column density, which they estimate to be N H ~ 3 x 10 24 cm 2 
with an uncertainty envelope extending a factor of three times lower 
and arbitrarily higher. The infrared luminosity derived on the basis 
of our SPIRE and PACS observations (see Figure|5]and Table[4} is 
2.0 x 10 13 L o . It appears to be fairly well isolated in the infrared, the 
nearest comparably bright 250 micron source being located half an 
arcminute away. 

Some of these ultraluminous galaxies could be lensed but a 
first look has not yet revealed these in our sample. The GOODS- 
N source (d), with coordinates (J2000.0) 189.309509 / 62.20232, 
is cited by Barger et al. (2008) as having a spectroscopic red- 
shift z = 3.157. Again, ambient nearby sources have relatively 
weak SPIRE fluxes unlikely to appreciably affect the SPIRE flux 
attributed to our source of primary interest. Another source only 
~ 10 arcsec away is also noted in the NASA/IPAC Extragalac- 
tic Database (NED). This appears not to have a measured 24 yum 
flux and is not listed by Barger et al. (2008 ). However, |Law et aT] 
(2007 ) have included this source in their discussion of distant ir- 
regular galaxies. The object designated as BX 150 appears elon- 
gated roughly along a north/south direction, is ~ 0.5 arcsec long 
with an aspect ratio roughly 2:1, and has a redshift z = 2.28. At 
optical wavelengths, it is 1.3 magnitudes fainter than the ultralumi- 
nous infrared source and, at its rather high displacement of ~ 10 
arcsec, it is unlikely to provide significant lensing. With this pro- 
viso, the intrinsic infrared luminosity of GOODS (d) appears to be 
~ 1.1 x 10 13 L o . 

Figure [5] provides the SEDs of these ultraluminous sources. 



6 D. Brisbin et al. 



Our combined surveys of GOODS-N and LN cover ~ 0.47 square 
degrees, or one part in 85,000 of the sky. Given that we observe sev- 
eral high-luminosity sources in the small area covered, it suggests 
that approximately 10 6 sources in the infrared luminosity range 
~ 10 13 L o should be observable in the Universe, at redshifts z = 2.5 
to 3 at the current epoch. The number of comparably luminous 
sources observable at lower redshifts appears to sharply decline. 

It is unlikely that three selection effects inherent in our obser- 
vations cast this conclusion into serious doubt: 

(i) The first is that Figure|4] on which the conclusion is based, 
only includes sources detected at all three SPIRE wavelengths. For 
highly redshifted sources, the 100 [im infrared emission peak is 
redshifted into the 500 yum range, favoring the detection of galaxies 
at all three wavelengths, including 500 pm. 

(ii) However, compensating for this effect, lower redshift 
sources are more readily detected by a factor inversely proportional 
to luminosity distance squared. Although these two effects partially 
cancel, high-luminosity sources should be more readily detected at 
low than at high redshifts. 

(iii) The XID catalogs search for SPIRE sources solely at loca- 
tions where Spitzer 24 yum sources exhibit flux densities 3> 20 /jJy in 
the GOODS-N field and > 50 yujy in LN. We may thus be missing 
sources at redshifts at which poorly emitting spectral regions are 
redshifted to 24/im. At redshifts z ~ 1.4, for example, the 9.7 /an 
silicate absorption dip shown by Spoon et al. ( 2007 1 to be prevalent 
in many ULIRGS is shifted to 24 //m. This may account for the 
striking absence of low-luminosity sources, at z ~ 1.4, i.e. the lack 
of sources hugging the luminosity distance curve at this redshift in 
Figure]?] 



7 DISCUSSION 

Far-infrared surveys with Herschel need to take source confusion 
and source blending into account, particularly at the longest wave- 
lengths, 350 and 500 yum. GOODS-N and LN are the two deepest 
surveys undertaken as part of the HerMES project to date. In these 
deep fields, crowding of sources presents especially serious prob- 
lems. In establishing a set of criteria that assess source blending, 
we have taken a preliminary step toward estimating the utility of 
the survey data for different purposes. This has proven itself use- 
ful in our analysis of the ultraluminous galaxies, some of which 
we described in Section 6 and whose characteristics are exhibited 
in Figure [5] and Table [4] In view of the high infrared luminosities 
we find, it is particularly satisfying that seven of the eleven sources 
cited turn out to be triply secure, i.e. with high purities in all three 
SPIRE wave-bands, and that five of the sources also are observed 
by PACS where blending is not severe, particularly in the 100 yum 
waveband. In this context, we have placed no upper limit on ac- 
ceptable values of II, which are especially high for GOODS(c) and 
LN(f), suggesting especially low ambient source contributions at 
their locations. 

In compiling the SEDs for GOODS-N and LN, we have 
elected to work with the S&K models because they are based on 
a limited set of well-defined physical parameters. The models thus 
make predictions that our SEDs may be able to verify, refute, or ex- 
tend. S&K do not specifically address the effects of adding an AGN 
component to a starburst model. However, they do provide a star- 
burst fit for NGC 6240 and propose that addition of a small AGN 
component could provide an improved fit. Most starbursts generally 
also exhibit some AGN activity. Perhaps because of this, the S&K 
models appear to provide reasonable fits. The major weakness of 



Table 5. Ratios (R = FfinlF opl j ca i) of SPIRE flux densities (consistently 
measured at 250 /im to minimize source blending) to optical flux densities 
at a rest wavelength ~3000 A, respectively measured at 4500 A, 6100 A, 
and 7600 A for successively larger redshifts. These three wavelengths were 
chosen to represent the visible spectra and avoid the Balmer jump at 3650 
A. For each redshift bin we compare the flux density ratio for the N 12 most 
luminous galaxies {RBright) to that of the N/2 least luminous galaxies (RDtm) 
where N is the total number of L > 10 12 L o sources. 



z 


N 


RBright 


RDim 


^Bright 1 Roim 


0.95 < z < 1.05 


18 


64500 


72000 


0.90 


1.9 < z < 2.1 


18 


62500 


17600 


3.6 


2.4 < z < 2.8 


10 


62300 


5920 


11 



the S&K models, as well as that of all others, tends to be the diffi- 
culty in accounting for the seemingly random relationship between 
the infrared and optical portions of the SEDs that is so apparent in 
Figures[2][3] and|5] 

We investigated the relationship between flux-density ratios at 
optical and far-infrared wavelengths in high- and low-luminosity 
galaxies. Current theory suggests that starbursts involve stellar 
mass distributions obeying the Salpeter initial mass function ( |Zin-| 
|necker & Yorke|2007| >. The drop in luminosities from the most mas- 
sive O type stars with mass ~ 120M o to the early B type stars at 10 
M Q , can then be shown to be roughly in a ratio of 500:1, i.e. with 
a contrast considerably higher than that of the mass ratio, roughly 
12:1. The highest mass ranges will thus be depleted most rapidly, 
ending their lives in supernova explosions in which at least some 
of the dust will be destroyed or expelled from the galaxy. The most 
luminous galaxies found using our SEDs and their associated red- 
shifts would thus be expected to be the very youngest as well as 
those most densely shrouded by dust, i.e. having the lowest frac- 
tional optical luminosities. To test this hypothesis we restricted our- 
selves to galaxies with total bolometric luminosities, L > 10 12 L o , 
as these have long been considered likely starburst mergers, albeit 
with potential contributions from AGNs ( Sanders & Mirabel 19961. 
In Table[5] we compare the flux density ratios for the most and least 
luminous sources in three redshift bins. It is evident that larger ra- 
tios correspond to more luminous sources at all redshifts but that 
these differences greatly diminish at lower redshifts. This finding is 
both new and significant. It indicates evolutionary trends that may 
need to be incorporated into more advanced models of starbursts 
designed to yield SEDs which not only mirror observed ratios of 
optical to infrared emission, but also define a galaxy's place in its 
evolutionary history. 

Figure [4] provides a capsule history of galaxy evolution over 
cosmological epochs for the sample included in our two deep sur- 
veys. The shapes of these distributions are nearly identical in LN 
and GOODS-N, which motivated us to plot both in the same figure. 
A glance confirms that sources having the highest luminosities are 
found at highest redshifts, i.e. earliest epochs. Luminosities higher 
than 10 13 L o are generally observed at redshifts z ~ 2 - 3.2, the 
highest redshifts reached in our surveys. 



8 CONCLUSIONS 

Confusion, which can be troubling at 250 yum, becomes increas- 
ingly severe at 350 and 500 yum. Yet the data at these longer wave- 
lengths are particularly important given how little is known about 
this spectral domain. We believe that the triply-secure sources listed 
in Tables 2, 3, and 4 will be in demand for follow-on studies that X- 
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ray astronomers, spectroscopists, and others may wish to undertake 
on sources known to be especially free of confusion. 

With SPIRE photometry data in hand along with cross identi- 
fications at several shorter wavelengths, we have constructed SEDs 
for a handful of trustworthy sources in the GOODS-N and LN re- 
gions. Many of these can be fit by starburst SED models, such as 
those created by Siebenmorgen & Kriigel (2007), which yield in- 
formation on luminosity, dust mass, and size. Figures 4 and 5 show 
a number of ultraluminous galaxies with L ;R ~ 10 13 L o . Although 
these are extreme systems, they do not appear to deviate from the 
general distribution at high redshift. A major strength of the deep 
HerMES surveys is their ability to obtain reliable source luminosi- 
ties and star-formation rates based on flux densities in the infrared 
and at auxiliary wavelengths as well as redshifts compiled in the 
XID catalogues. 
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Figure 3. GOODS-N galaxy SEDs plotted as in Figure [I] The 500 iim 
measurements for GOODS sources A, C, E, F, H, K, and L are compatible 
with zero, so here we plot only the upper error bar, barely visible within the 
diamond symbol. As in Figure|2] we weighted the 24 (im data a quarter as 
heavily as the longer wavelengths. Observational data shortward of 24 /jm 
is plotted for reference but not used in fitting. 
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Infrared Source Luminosities 




redshift (z) 



Figure 4. Infrared source luminosities (integrated over 8-1000 fim) in 
GOODS-N and LN plotted as a function of redshift for all sources detected 
at all three SPIRE wavelengths. Diamonds indicate GOODS-N luminosities 
obtained from SEDs fitted by S&K models, and triangles indicate similar 
luminosities for LN sources. The solid line shows the growth of luminosity 
distance squared with z. It serves as a rough lower bound to the luminosi- 
ties in our selection of observed sources; the scatter of data points about the 
curve provides a visual impression of the uncertainties in those luminosities. 
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Figure 5. Spectral energy distributions of the most luminous sources in 
GOODS North (a) to (d) Lockman North (e) to (k). As in Figures 2 and 3, 
we solely plot the upper error bars at 500 /im for sources (a) to (e), whose 
lower error bars are compatible with zero. We similarly plot an upper error 
bar at 170 /jm for LN (k). Luminosities and positions for these sources are 
presented in Table|4] Seven of these eleven sources are triply secure; two of 
these are also included in Figures 2 and 3. A brief description of sources (c), 
(d), (f) and (i) is provided in Section 6. Although some of the ultraluminous 
sources exhibit significant AGN activity, we have applied S&K model fits, 
as discussed in Section 7. As in Figures 2 and 3 we have weighted the 24 /im 
data only one quarter as heavily as the longer wavelengths. Observational 
data shortward of 24 (im is plotted for reference but not used in fitting. 
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Table 2. Data on LN Galaxies whose SEDs appear in Figure [2] We list the ID as established by Roseboom et al. ( 201 0;t, right ascension and declination 
(J2000), photometric redshift z > 0.5, our purity indices LL;, total bolometric luminosity as estimated by an S&K model, star formation rate based on the 
infrared luminosity relation ( Kennicutt 1998 1, and dust mass estimated by the S&K model. Along with these derived parameters, we list the model parameters: 
nuclear radius, visual extinction to center (A,,), and gas density within hotspots (m, s ). 



ID 


RA 


Dec 




^-boiometric C^o) 


SFR (M Q /yr) 


M dust ( M o) 


n 250 


n 350 


n 500 


Radius (kpc) 


Ay 


L °gl0H.r xcm3 ) 




LN (A) 


161.8365 


59.1211 


0.56 


4.0X10 11 


66 


1.6xl0 8 


0.89 


0.67 


0.74 


3.0 


35.9 


3 


LN (B) 


161.5001 


58.8732 


0.60 


7.9X10 11 


130 


3.2xl0 8 


0.74 


0.62 


0.88 


3.0 


72.0 


4 


LN(C) 


161.1277 


59.1956 


0.72 


1.3X10 12 


210 


3.2x10 s 


1.10 


1.00 


1.07 


3.0 


72.0 




LN (D) 


161.3232 


59.2086 


0.82 


I.3X10 12 


210 


3.2x10 s 


0.88 


0.80 


0.98 


3.0 


72.0 


4 


LN (E) 


161.0530 


59.0762 


0.94 


2.0X10 12 


330 


1.6x10 s 


0.92 


0.65 


0.69 


3.0 


35.9 


3 


LN (F) 


161.3680 


59.2242 


0.96 


2.5X10 12 


400 


3.6x10 s 


0.96 


0.97 


0.56 


9.0 


9.0 


4 


LN (G) 


161.3429 


59.2269 


1.02 


4.0xl0 12 


650 


1.6x10 s 


0.89 


0.87 


0.64 


3.0 


35.9 


2 


LN (H) 


161.4871 


58.8886 


1.06 


2.5xl0 12 


410 


1.6x10 s 


0.76 


0.82 


0.93 


3.0 


35.9 


2 


LN 0) 


161.8669 


58.8708 


2.76 


2.5xI0 13 


4300 


4.9x1 9 


0.91 


0.96 


0.77 


9.0 


120,0 


4 



Table 3. Data on GOODS-N galaxies whose SEDs appear in Figure[5] listed as in Table[2] The redshifts here are spectroscopic. 



ID 


RA 


Dec 




L bolometric (^o) 


SFR (M /yr) 


M dust 


n 250 


n 350 


n 500 


Radius; (kpc) 


A v 


L °gl0< n /i J x 




GOODS (A) 


189.0274 


62.1643 


0.6380 


5.0X10 11 


84 


3.2X10 8 


0,83 


0,84 


0.44 


3.0 


72.0 


4 




GOODS (B) 


189.3938 


62.2898 


0.6402 


2.0xl0 12 


340 


6.0x1 7 


0,98 


0,64 


0,90 


1.0 


119.0 


2 




GOODS (C) 


189.2979 


62.1820 


0.8549 


1.3xl0 U 


21 


3.2x10 s 


1.00 


0,81 


0,74 


3.0 


72.0 


4 




GOODS (D) 


189.1403 


62.1683 


1.0160 


1.6xl0 12 


260 


1.6x10 s 


0.99 


0,75 


0.54 


3.0 


35.9 


3 




GOODS (E) 


189.0633 


62.1691 


1.0270 


1.3xl0 12 


210 


1.6x10 s 


0.92 


0,52 


0,45 


3.0 


35.9 


2 




GOODS (F) 


189.3171 


62.3541 


1.1440 


l.OxlO 12 


160 


S.lxlO 7 


0,73 


0.68 


0.68 


3.0 


17.9 


2 




GOODS (G) 


189.1438 


62.2114 


1.2242 


2.5xl0 13 


4300 


6.0xl0 7 


0.88 


0.89 


0.98 


1.0 


119.0 


4 




GOODS (H) 


189.2137 


62.1810 


1.2258 


6.3X10 11 


100 


1.6xl0 8 


0.71 


0.83 


0.31 


3.0 


35.9 


4 




GOODS (I) 


189.2614 


62.2338 


1.2480 


1.3xl0 12 


200 


3.6xl0 8 


0.98 


0.93 


0.97 


9.0 


9.0 


4 




GOODS (J) 


189.2566 


62.1962 


1.7600 


6.3xl0 12 


1000 


7.3x10 s 


1.02 


1.04 


1.03 


9.0 


18.0 


4 




GOODS (K) 


189.3036 


62.1955 


1.8150 


2.0xl0 12 


310 


S.lxlO 7 


1.08 


0.86 


0.50 


3.0 


17.9 


3.4 




GOODS (L) 


189.0764 


62.2640 


2.0000 


2.0xl0 13 


3400 


6.0x1 7 


0.82 


0,71 


0.33 


1.0 


119.0 


2 





Table 4. Data on luminous galaxies whose SEDs appear in Figure[5] listed as in previous tables. The GOODS sources have spectroscopic redshifts, as do LN 
(f) and (j). The rest of the LN sources have photometric redshifts. The sources GOODS (a) and LN (i) correspond to the sources GOODS (L) and LN (I) shown 
in the previous figures and tables. 



ID 


RA 


Dec 




L bolometric t 1 ©) 


SFR (M /yr) 




n 250 


n 350 


n 500 


Radius (kpc) 


A v 


L °gl0( n fti x 


m 3 ) 




GOODS (a) 


189.0764 


62.2640 


2.0000 


2.0xl0 13 


3400 


6.0x1 7 


0.82 


0.71 


0.33 


1.0 


119.0 


2.0 




GOODS (b) 


189.1483 


62.2400 


2.0050 


2.5X10 13 


4200 


3.6xl0 7 


0.69 


0.24 


0.38 


1.0 


71.0 


2.0 




GOODS (c) 


188.9901 


62.1734 


3.0750 


6.3X10 12 


1000 


2.0x1 9 


1.61 


1.22 


0.84 


15.0 


18.0 


4.0 




GOODS <d) 


189.3096 


62.2024 


3.1569 


1.3X10 13 


2000 


4.0x1 7 


0.57 


0.55 


0.84 


3.0 


9.0 


2.0 




LN (e) 


161.7059 


59.3247 


1.28 


I.3XI0 13 


2100 


3.2x10 s 


0.99 


0.85 


0.20 


3.0 


72.0 


2.0 




LN (f) 


161.0415 


58.8735 


2.28 


2.5xl0 13 


3400 


4.1xl0 5 


0.98 


1.10 


1.71 


0.3 


6.7 


2.0 




LN <g) 


161.5408 


58.7950 


2.58 


I.3XI0 13 


2000 


2.0x1 9 


0.82 


0.95 


0.58 


15.0 


18.0 


4.0 




LN (h) 


160.9635 


58.9555 


2.68 


2.5xl0 13 


4300 


8.8xl0 6 


0.18 


0.95 


0.96 


0.3 


144.0 


2.0 




LN (i) 


161.8667 


58.8704 


2.76 


2.5xl0 13 


4300 


4.9xl0 9 


0.91 


0.96 


0.77 


9.0 


120.0 


4.0 




LN (j) 


161.5541 


58.7886 


2.96 


2.0xl0 13 


3100 


4.0x1 7 


0.76 


0.61 


0.60 


3.0 


9.0 


3.0 




LN (k) 


161.8259 


59.2771 


3.12 


I.3XI0 13 


2000 


8.0xl0 7 


0.98 


0.70 


0.54 


3.0 


17.8 


2.0 





